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Abstract 

The generalized lambda distribution (GLD) is a flexible four parameter distribution 
with many practical applications. L-moments of the GLD can be expressed in closed 
form and are good alternatives for the central moments. The L-moments of the GLD 
up to an arbitrary order are presented, and a study of L-skewness and L-kurtosis 
that can be achieved by the GLD is provided. The boundaries of L-skewness and L- 
kurtosis are derived analytically for the symmetric GLD and calculated numerically 
for the GLD in general. Additionally, the contours of L-skewness and L-kurtosis are 
presented as functions of the GLD parameters. It is found that with an exception 
of the smallest values of L-kurtosis, the GLD covers all possible pairs of L-skewness 
and L-kurtosis and often there are two or more distributions that share the same 
L-skewness and the same L-kurtosis. Examples that demonstrate situations where 
there are four GLD members with the same L-skewness and the same L-kurtosis are 
presented. The estimation of the GLD parameters is studied in a simulation example 
where method of L-moments compares favorably to more complicated estimation 
methods. The results increase the knowledge on the distributions that belong to the 
GLD family and can be utilized in model selection and estimation. 
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1 Introduction 



The generalized lambda distribution (GLD) is a four parameter distribu- 
tion that has been applied to various problems where a flexible parametric 
model for univariate data is needed. The GLD provides fit for a large range 
of skewness and kurtosis values and can approximate many commonly used 
distributions such as normal, expon ential and unif orm. The applications of 
the GLD include e.g. option pri cing (jCorradd . l200ll ). indepe ndent com ponent 
analysis (IKarvanen et all I2002T). statistical proces s control (IPall . 120051 ) . anal- 



ysis of fatigue of ma terials (IBigerelle et al.l . 120051 ) . meas urement technology 



( Lampasi et all 120051) and generation of random variab les (IRamberg and Schmeiserl . 
1974J : lKarvapfi3 . 12OO3I : iHeadrick and MugdadibL l2006h . 



The price for the high flexibility is high complexity. Probability density func- 
tion (pdf) or cumulative distribution function (cdf) of the GLD do not exist in 
closed form but the distribution is defined by the inverse distribution function 
( IRamberg and Schmeiserl . Il974l ) 
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where < u < 1 and Aj, A2, A3 and A4 are t he parameters of th e GLD. 
Equation (1) defines a distribution if and only if (IKarian et al.l . Il996l ) 
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— > for all u G [0,1]. 



(2) 



Karian et al.l (119961 ) divide the GLD into six regions on the basis of the param- 
eters A 3 and A 4 that control skewness and kurtosis of the distribution. These 
regions are presented in Figure 1. The regions have different characteristics: 
the distributions in region 3 are bounded whereas the distributions in region 4 
are unbounded; the distributions in regions 1 and 5 are bounded on the right 
and the distributions in regions 2 and 6 are bounded o n the left. The bound- 
aries o f the domain in each regio n are described e.g. by lKarian and Dudewicz 
(l2000h and lFournier et all (l2007h . 



The members of the GLD family have been traditionally characterized by 
the central moments, especially by the central momen t skewness a nd kurto- 



sis (|Ramberg and Schmeiserl . Il974j ; IKarian et allll996l ). L-moments (IHoskingl . 
19901 ) . defined as linear combinations of order statistics, are attractive al- 



ternatives for the central moments. Differently from the central moments, 
the L-moments of the GLD can be expressed in closed form, which allows 
us to derive some analytical results and makes it straightforward to per- 
form numerical analysi s. The parameters of t he GLD can b e estimated by 
method of L-moments (IKarvanen et al.l . 120021 ; lAsquithl . 120071 ). Method of L- 
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A 4 
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3 - 2 
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Fig. 1. GLD regions as denned in dKarian et all llflfld ). The distributions are 



bounded on the right in regions 1 and 5, bounded on the left in regions 2 and 
6, bounded in region 3 and unbounded in region 4. 



moments can be used independently or together wi th other estimation m eth- 



ods, analogously to wh a t was done in (jSul . 120071 ) , (IFournier et all 120071 ) and 



( Lakhany and Mausser . 2000l ) with numerical maximum li k elihoo d, method 
of moments, m ethod of percentiles ([Karian and Dudewiczl . Il999l). the least 
squar e method (jOztiirk and Dald . ll985l ) and the starship method (IKing and MacGillivray 
1999h . 



Karian and Dudewiczl (120031 ) report that method of percentiles gives superior 
fits compared to method of moments. They also study method of L-moments 
in some of their examples where the overall performance looks comparable to 
the overall performance of the percentile method. Method of percentiles and 
method of L-moments are related in sense that they both are based on order 
statistics. There are, however, some differences: L-moments are linear combina- 
tions of all order statistics whereas method of percentiles uses only a limited 



numb er of order statistics that need to be explicitly chosen. IFournier et ah 



(120071 ) study the choice of the order statistics for the percentile method and 
conclude that there is no trivial rule for choosing them. Such problems are 
avoided in method of L-moments. 



In this paper we analyze the GLD using L-moments and consider the estima- 
tion by method of L-moments. In Section 2 we present the L-moments of the 
GLD up to an arbitrary order. In Section 3 we analyze the symmetric case 
(A 3 = A 4 ) and derive the boundaries of L-kurtosis analytically and in Section 4 
we consider the general case (A3 7^ A4) and calculate the boundaries of the 
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GLD in the terms of L-skewness and L-kurtosis using numerical methods. The 
boundaries are calculated separately for each GLD region. We also calculate 
the contours of L-skewness and L-kurtosis as functions of A3 and A4. The ex- 
amples presented illustrate the multiplicity of the distributional forms of the 
GLD. In Section 5 we present an example on the est imation by me t hod o f 



L-moments and compare the results to the results by iFournier et al.l (120071 ). 
Section 6 concludes the paper. 



2 L-moments and the GLD 



The L-moment of order r can be expressed as ( IHoskingi . Il990h 

JO 

where 

Ec-ir*- 1 



p:-i{u) 



r — l\(r + k — \\ k 



k=0 



(3) 



(4) 



is the shifted Legendre polynomial of order r — 1. All L-moments of a real- 
valued random variable exists if and only if the random variable has a finite 
mean and furthermore, a distribu t ion whose m ean exists, is uniquely deter- 
mined by its L-moments (IHoskingi . Il990l . 120061 ) . 



Similarly to the central moments, L% measures location and L 2 measures scale. 
The higher order L-moments are usually transformed to L-moment ratios 



r = 3,4, 



(5) 



L-skewness T3 is related to the asymmetry of the distribution and L-kurtosis 
T4 is related to the peakedness of the distribution. Differently from the central 
moment s kewness and kurto sis, r 3 and r 4 are constrained by the conditions 
flHoskind . Il990l : Ijonesl . l2004h 



-1 <r 3 < 1 



(6) 



and 



(5r 3 2 - l)/4 <r 4 < 1. 



(7) 



Distributions are commonly characterized using an L-moment ratio diagram 
in which L-skewness is on the horizontal axis and L-kurtosis is on the vertical 



axis. 
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The L-moments of the GLD have been presented in the lite rature up to or- 
der 5 (IBergevinl . Il993l ; iKarvanen et all 120021 ; lAsquithl . 120071 ). We generalize 
the results for an arbitrary order r. The derivation is straightforward if we 
first note that 
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where the last equality follows from the property P*_ 1 (l— u) 
We obtain for r > 2 



\ 2 L r = J2(-iy 

k=0 
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(9) 



k+l + \ 3 k + 1 + A 4 



The explicit formulas for the first six L-moments of the GLD are 



1 



Ll - Al "A^Vl + A 4 



L0X0 — — 
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L,A 
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I + A3 2 + A 3 3 + A 3 I + A4 2 + A4 3 + A 4 
1 12 30 20 

La\ 2 = ~ ; , ; + , ; + 



I + A3 
1 



2 + A 3 3 + A 3 
12 30 



4 + A 3 
20 



5^2 



I + A4 2 + A4 3 + A4 4 + A4 
1 20 90 140 70 
+ t: ; z r- + 



1 + A 3 2 + A 3 3 + A 3 

1 20 90 

+ 



4 + A 3 ' 5 + A 3 
140 70 



I + A4 2 + A4 3 + A4 4 + A4 



5 + A 4 



(10) 



(11) 
(12) 

(13) 
(14) 
(15) 



and 



1 



6^2 



+ 



30 



210 560 

+ 



630 252 
+ 



I + A3 2 + A 3 3 + A 3 4 + A 3 5 + A 3 6 + A 3 
1 30 210 560 630 252 

1 + A 4 2 + A 4 3 + A 4 4 + A 4 5 + A 4 6 + A 4 ' 



(Asquith, 


2007 


) and ( 


Bergevin , 


1993) 



(16) 



but expressed in a different form. 
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The mean of GLD and therefore also all L- moments exist if A3, A4 > —1. This 
implies that the characterization by L-moments covers regions 3, 5 and 6 and 
the subset of region 4 where A 3 , A 4 > —1. 



3 The L-moments of the GLD in the symmetric case 



In this section we consider the special case A3 = A4, which defines a symmetric 
distribution. The symmetric distributions can be bounded (in region 3) or 
unbounded (in region 4). Applying the condition A3 = A 4 to equations (13) 
and (14) we find that T3 = and 



4 



3A3 + 2 



A| + 7A 3 + 12 : 



A3 



A 4 > -1. 



(17) 



In Figure 2 r 4 is plotted as a function of A 3 . The value of r 4 approaches to 1 
as A3 approaches to infinity. The minimum of r 4 is found to be 



12-5^6 
12 + 5a/6 

and is obtained when A3 = — 1 + \^6. 



-0.0102, 




A3 

Fig. 2. L-kurtosis t^A^) for symmetric distributions A3 



A 4 . 



In the symmetric case it is also possible to solve A 3 and A 4 as functions of r 4 
(IKarvanen et al.l . |2002| ) 



A 4 — A3 



3 + 7r 4 ± y/1 + 98r 4 + ri 
2(1 - r 4 ) 



(19) 



It follows that if T4 > 1/6 there are symmetric GLD members available from 
both regions 3 and 4. If 

12-5^ 1 , N 

< n < -, (20) 



12 + 5\/6 



6' 
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there are two GLD members available from region 3, and if 



r 4 < 



12-5^6 
12 + 5^' 



(21) 



there are no symmetric GLD members available. 

Examples of symmetric GLDs sharing the same r 4 are presented in Figure 3. 
In Figure 3(a) the value of r 4 is set to be equal to the r 4 of the normal 
distribution. Both GLDs are from region 3 and are bounded. The GLD for the 
more peaked distribution in Figure 3(a) has r 6 m 0.0004 whereas the other 
GLD has r@ w 0.043, which is very near to the Tq of the normal distribution. 
For comparison the pdf of the normal distribution is plotted (dotted line). The 
differences with the normal distribution are appreciable only in the far tails. 
In Figure 3(b) r 4 = 0.25, which means that equation (17) has one positive 
and one negative root. The GLD with tq « 0.016 and the higher peak is from 
region 3 and has a bounded domain. The other GLD with tq « 0.121 is from 
region 4 and has unbounded domain. 
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Fig. 3. Examples of symmetric GLDs sharing the same r 4 . The pdfs of the GLDs 
are plotted with solid line and the pdf of normal distribution is plotted with dotted 
line. All distributions have Ai = 0, A2 = 1 and T3 = 0. 



In addition to the symmetric case we may consider some other special cases. 
If A3 = we obtain 

r 3 (X t ) =^ (22) 
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and 



If A a = we obtain 



r s (A 3 ) =^2— ^ (24) 

A 3 + <J 

and 

Interestingly, equations (23) and (25) are identical to equation (17), which 
indicates that the (A 3 , A 4 ) pairs (A, A), (A, 0) and (0, A) lead to the same value 
of r 4 . 



4 Boundaries of the L-moment ratios of the GLD 



In the general case it is difficult to derive analytical results and therefore we 
resort to numerical methods. Because the L-moments of the GLD are available 
in closed form, we may choose a straightforward approach and calculate the 
values of r 3 and r 4 for a large number of (A 3 , A 4 ) pairs. Naturally, the values 
of (A 3 , A 4 ) need to be chosen such a way that the essential properties of the 
GLD are revealed. In region 4, we use a grid of one million points where both 
A3 and A 4 are equally spaced in the interval [—1 + 1CT 10 , lCT 10 ]. In region 3 
we use an unequally spaced grid of 4104676 points. The grid is dense near 
the origin (A3 and A 4 have minimum 1CT 10 ) and sparse for large values (A 3 
and A 4 have maximum 10 12 ). In region 5, A3 is equally spaced in the interval 
[— 1 + 10~ 10 , 0] and A 4 is unequally spaced starting from —1. The grid contains 
1362429 valid points. The calculations thus result in a large number of (r 3 , r 4 ) 
pairs. The results for region 6 can be obtained from the results for region 5 
by swapping A3 and A 4 . The problem is to find the boundaries of the GLD in 
the (t3,t 4 ) space on the basis of these data. Let (r 3 [i], T^[j\) be the image of 
the grid point (A 3 [i], A 4 [j]), where i and j are the grid indexes. If the point 
(r 3 [z], r 4 [j]) is not located inside the polygon defined by points (r 3 [z — 1], r 4 [j]), 
+ 1]), (r 3 [i + l),n\j)) and (r 3 [i],r 4 [j - 1]) it is considered to be a 
potential boundary point. The images of the points that are located on the 
boundary of the grid in the (A 3 , A 4 ) space are also potential boundary points 
in the (r 3 , r 4 ) space. The actual boundaries are found combining the potential 
boundaries and finally it is checked that all (r 3 , r 4 ) pairs are located inside the 
boundaries. The calculations are made using R (|R Development Core Team! 
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20061) and the co ntributed R packages PBSmapping (ISchnute et all 120061 ) and 
gld ( IKingl . 120061 ) are utilized. 



The boundaries are shown in Figure 4. It can be seen that region 3 covers most 
of the (r 3 , r 4 ) space; only the smallest values of r 4 are unattainable. Region 4 
largely overlaps with region 3 but cannot achieve r 4 smaller than 1/6. Region 
5 and its counterpart, region 6, cover a rather small area of the (r 3 , r 4 ) space. 




-1.0 -0.5 0.0 0.5 1.0 -1.0 -0.5 0.0 0.5 1.0 

(c) Region 5 (d) Region 6 

Fig. 4. Boundaries of the GLD in the (r 3 , r 4 ) space. The shading shows the values of 
L-skewness and L-kurtosis that are achievable by the GLD. The outer limits present 
the boundaries for all distributions. 

Figure 5 shows the (r 3 , r 4 ) area where there exist GLD members from regions 3, 
4 and 5, or from regions 3, 4 and 6. To illustrate these distributions we fix Ai = 
0, A2 = 1, t 3 = 0.4 and r 4 = 0.25 and seek for the solutions in regions 3, 4 and 
6. Table 1 and Figure 6 present the four distributions that are found. Although 
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Table 1 

Four different GLD members with L\ = 0, L 2 = 1, T3 = 0.4 and T4 = 0.25. 



Region 


Ai 


A 2 


A3 


A 4 


T5 


T6 


3(a) 


5.322 


0.138 


21.526 


0.286 


0.163 


0.103 


3(b) 


-1.168 


0.124 


5.417 


92.608 


-0.029 


0.067 


4 


-1.62 


-0.157 


-0.014 


-0.212 


0.158 


0.121 


6 


-7.04 


-0.194 


11.905 


-0.306 


0.204 


0.180 



the first four L-moments of the distributions are the same, the differences of 
the pdfs are clearly visible. Distribution 3(a) has bounded domain with sharp 
left limit. Distribution 3(b) is also bounded but characterized by the high 
peak. Distribution 4 has unbounded domain but is otherwise almost similar 
to distribution 3(a). Distribution 6 is bounded from the left and characterized 
by a minor 'peak' around x — 2. 
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-1.0 -0.5 0.0 0.5 1.0 

x 3 

Fig. 5. The (t 3 ,T4) area where there exist GLD members from regions 3, 4 and 5 
(negative L-skewness), or from regions 3, 4 and 6 (positive L-skewness). 

The contours of T3 and T4 in different GLD regions are shown in Figures 7, 8 
and 9. In region 3 the contours are rather complicated. From Figure 7(a) and 
7(c) it can be seen that there are three separate areas where r 3 has negative 
values. It is also seen that besides the line A3 = A4, there are two curves where 
T3 = 0. Naturally, the distributions defined by the curves are not symmetric; 
they just have zero L-skewness. 



5 Estimation by method of L-moments 

In method of L-moments, we first calculate L-moments L 1 , L 2 , f 3 and f 4 
from the data. Then we numerically find parameters A3 and A 4 that minimize 
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Fig. 6. Pdfs of four GLDs with t% = 0.4 and T4 = 0.25. The parameters of the 
distributions are presented in Table 1. 

some objective function that measures the distance between the L-skewness 
and L-kurtosis of the data (73,74) and the L-skewness and L-kurtosis of the 
estimated model (73 (A3, A 4 ), 7 4 (A 3 , A 4 )). After that estimates A\ and A 2 can be 
solved from equations (11) and (12). A natural objective function is the sum 
of squared distances 



(73 -7 3 (A 3 ,A 4 )) + (74 - 7 4 (A 3 , A 4 )) 



(26) 



which was used also by lAsquithl (120071 ) 



We st udy the same simulation example that was studied by iFournier et al 



(120071 ). They compare five estimators in situation where a sample of 1000 
observations is generated from GLD(0, 0.19, 0.14, 0.14), which is a symmet- 
ric distribution close to the standard normal distribution. We apply method 
of L-moments to the same problem and present the results in a compara- 
ble form. Before carrying out the simulations, we analyze the example us- 
ing the results given in the preceding sections. Because A 3 = A 4 = 0.14, 
GLD(0, 0.19, 0.14, 0.14) is a symmetric distribution from region 3. The the- 
oretical L-moments L x = 0, L 2 « 0.60407, 7 3 = and 7 4 « 0.12305 of 
GLD(0, 0.19, 0.14, 0.14) are obtained from equations (11)-(14). Because 7 3 = 
we may apply equation (19) that gives the "correct" solution A3 = A4 = 0.14 
and another solution A3 = A4 ~ 4.26316 also from region 3. By inspect- 
ing the contours in Figure 7 we find that there exist also two asymmetric 
GLDs from region 3 with 7 3 = and 7 4 « 0.12305. The numerical mini- 
mization of objective function (26) with suitable initial values gives solutions 
(A 3 « 1.98, A 4 « 22.59) and (A3 « 22.59, A 4 « 1.98). Figure 10 shows the 
contours of objective function (26) when 73 = and 74 ~ 0.12305. It can be 
clearly seen that the achieved numerical solution depends on the initial values 
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(c) Contours of L-skewness T3 (zoomed) (d) Contours of L-kurtosis T4 (zoomed) 
Fig. 7. Contours of T3 and T4 in region 3. The shading indicates negative values of 

of A3 and A4. 

In the simulation we choose a strategy where f 4 computed from the data 
and the two solutions of equation (19) are tried as initial values of A3 and A4. 
The numerical optimi zation of the objective fu nctio n (26) is carried out by the 
Nelder-Mead method ( Nelder and Meadl . 1965 ) in R ( R Development Core Team . 
20061 ) a nd the computer used fo r the simulation is comparable to the computer 
used by iFournier et al.l (120071 ). The results from 10000 data sets of 1000 ob- 
servations generated from GLD(0, 0.19, 0.14, 0.14) are presented in Table 2. 
The first column of the table reports the L-moment estimates corresponding 
to the smaller initial solution which leads to the estimates c orresponding to 



the generating GLD. The last three columns are copied from (IFournier et al. 
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Fig. 8. Contours of T3 and T4 in region 4. The shading indicates negative values of 




2007I ) and are based on 100 data sets of 1000 observations generated from the 



same GLD. The L-moment estimates of Ai, A3 and A4 are unbiased and A2 is 
slightly underestimated. Method of L-moments is about 100 times faster than 
the Fournier et al. method and over 7000 times faster than the starship meth- 
ods. In terms of the Kolmogorov-Smirnov statistics Ejcs between the observed 
data and the fitted model, method of L-moments gave slightly better fit than 
the starship methods with the used settings. 



13 



i 1 1 1 1 r ~i i i i i r 

5 10 15 20 25 0.0 0.2 0.4 0.6 0.8 1.0 

(a) Contours and the four solutions (b) Contours in the neighborhood of 

GLD(0, 0.19, 0.14, 0.14) 

Fig. 10. Contours of objective function (26) when T3 = and T4 ~ 0.12305. The 
dots show the two symmetric and the two asymmetric solutions. 

6 Conclusion 



We have presented the L-moments of the GLD up to an arbitrary order and 
studied which values of L-skewness r 3 and L-kurtosis r 4 can be achieved by 
the GLD. For the symmetric case the boundaries were derived analytically 
and in the general case numerical methods were used. It was found that with 
an exception of the smallest values of T4, the GLD covers all possible (V3, T4) 
pairs and often there are two or more distributions sharing the same T3 and 
T4. The example in Section 4 demonstrates a situation where there are four 
GLD members sharing the same r 3 and T4. 



We argue that L-moments are natural descriptive statistics for the GLD be- 
cause they, differently from the central moments, can be expressed in closed 
form. The existence of the closed form presentation of the GLD L-moments 
follows from the fact that the inverse distribution function of the GLD is avail- 
able in closed form and can be integrated. The relation between L-moments 
and distributions d efined by the inve rse distribution function is not restricted 
only to the GLD (jKarvanenl . 120061 ). The existence of the first four central 
moments of the GLD requires that A3, A4 > —1/4 whereas the existence of 
L-moments of any order requires only that A3, A4 > —1. Thus, using the L- 
moments we can charact erize wider subset of the GLD than using the central 
moments. lAsquithl (120071 ) concluded that even wi der subset of the GLD cou ld 
be characterized using the trimmed L-moments (jElamir and Seheultl . 120031 ). 



The results presented in this paper can be utilized in model selection and 
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Table 2 

Comparison of estimation methods. The mean and the estimated standard error 
of the L-moment estimators are computed from 1000 data sets generated from 
GLD(0, 0.19, 0.14, 0.14). The results for the other methods are based 100 gener - 
ated data sets from the same GLD and are copied from ( Fournier et al. . 20071 ) . 
Fournier et al. method is based on the percentile method and the minimization of 
the Kolmogorov-Smirnov distance in the (A,s,A4 ) space. Starship KS and star ship 
AD methods correspond to the starship method ( King and MacGillivray . 19991 ) us- 
ing the Kolmogorov-Smirnov or the Anderson-Darling distance, respectively. 



Quantity 



L-moments Fournier et al. Starship KS Starship AD 



Ai 


Mean 


0.00153 


-0.00030 


0.25010 


0.00720 




Std error 


0.10287 


0.05010 


0.10000 


0.07300 


A 2 


Mean 


0.18795 


0.19960 


0.20020 


0.19400 




Std error 


0.03625 


0.05650 


0.04760 


0.03020 


A 3 


Mean 


0.14012 


0.15110 


0.15530 


0.14600 




Std error 


0.03602 


0.05150 


0.04630 


0.03170 


A 4 


Mean 


0.13981 


0.15060 


0.14750 


0.14350 




Std error 


0.03573 


0.05030 


0.04510 


0.02860 


Time (s) 


Mean 


0.02767 


3.19010 


206.99000 


213.47000 




Std error 


0.00575 


0.15020 


40.42010 


42.15210 


Ers 


Mean 


0.01597 


0.01850 


0.01680 


0.01856 




Std error 


0.00319 


0.00430 


0.00270 


0.00360 



estimation. The characterization by L-moments gives an insight into which T3 
and T4 are available in each GLD region. For instance, there are symmetric 
GLD members available from both region 3 and 4, if T4 > 1/6. This kind of 
information is useful when making the decision whether the GLD is a potential 
model for certain data. The results are also useful in the estimation of the 
GLD parameters. The parameters can be estimated directly by method of L- 
methods or the L-moment estimates can be used as starting values for other 
estimation methods. In both cases, we can use the L-moment ratio boundaries 
to specify the potential GLD regions where we should search for the parameter 
estimates. The choice between alternative solutions can be based on the type 
of the domain (bounded/unbounded) or some other additional criterion such 
as the Kolmogorov-Smirnov statistic or L-moment ratios r 5 and r@. 

In the simulation example, method of L-moments compared favorably to more 
complicated estimation methods. Estimation by method of L-moments was the 
fastest in the comparison and the estimates were unbiased or nearly unbiased. 
The goodness of fit measured by the Kolmogorov-Smirnov statistic was the 
same or better than with the alternative estimation methods. More simulations 
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are needed to find out how general these results are. 
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